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1 .  INTRODUCTION 


The  Nuclear  Damage  Assessment  Computer  Code  (NUDACC)^  described  in 
this  report  is  an  outgrotrth  of  a  requirement  of  the  U«S«  Army's  Theater 
Nuclear  Force  Survivability  Program  to  provide  accurate  assessment  of 
the  survival  of  blue  theater  nuclear  forces  after  a  USSR-Warsaw  Pact 
(red)  nuclear  attack.  Assessments  of  this  nature  have  often  been  made 
by  using  a  cookie-cutter  methodology.  A  circle  of  appropriate  radius  is 
drawn  furound  the  actual  ground  zero  (A6Z),  and  items  within  the  circle 
are  considered  killed. 

A  moment's  reflection  on  this  approach  will  confirm  that  although 
the  simplicity  of  method  is  attractive,  it  is  certainly  not  realistic. 

The  NUDACC  methodology  determines  a  probed>ility  of  survival  for 
personnel  and  equipment  based  on  a  cumulative  log-normal  function  of  a 
particular  nuclear  weapons  effects  (NWE)  environment.  M2my  equipment 
items  are  vulnerable  to  several  NWE  environments,  in  which  case  the 
Individual  probabilities  are  multiplied  to  determine  the  probability  of 
survival  of  the  item. 

To  simplify  further  discussion,  it  is  necessary  to  define  the 
following  terms: 

UNIT  Any  combination  of  personnel  and  equipment  that 

constitutes  a  logical  functional  entity  and  can  be 
contained  in  a  rectangular  area  as  small  as  50  m  on  a 
side  or  as  large  as  500  m  on  a  side  (This  size 
restriction  is  due  to  the  current  data  dimensions  of 
the  program  and  can  easily  be  modified.) 

UNIT  NAME  Name  assigned  to  unit 

UNIT  IDENTI-  Unique  seven-digit  integer  number  assigned  to  each  unit 

FICATION  (This  number  is  used  to  correlate  corresponding  data  in 

NUMBER  several  files.) 

UNIT  IfNGTH,  Geometric  parameters  of  unit  (in  meters):  LENGTH  =  dis- 
WIDTH  tance  from  front  to  rear;  WIDTH  =  distance  across  front 

(Each  value  C2ui  range  from  50  to  500  m.) 

RIGHT  REAR  Corner  to  right  and  rear  of  observer  located  in  cen- 

CORNER  ter  of  unit  €md  looking  toward  front  (This  corner  is 

used  to  locate  unit.) 

^J.  Michalowicz,  R.  Moore,  and  K,  Siteasy,  NUDACC — A  Nuclear  Damage 
Assessment  Computer  Code  (U),  Harry  Diamond  Laboratories  HDL-PR-78-3 
(November  1978),  (CONFIDENTIAL) 
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UNIT  LOCA¬ 
TION 


nie  x-y  location  of  unit's  right  rear  corner  (in  meters) 
(Origin  of  coordinate  system  must  be  chosen  so  that  all 
x-y  values  are  positive.) 


UNIT  ROTA-  Angle  in  degrees  about  right  rear  corner  measured  coun- 

TION  ANGLE  terclockwise  from  horizontal 


EQUIPMENT  List  of  personnel  and  equipment  organic  to  unit 

MENU 


EQUIPMENT  unique  three-digit  integer  number  used  to  identify  item 

CODE  in  equipment  menu 


DOMINANT  KILL  Most  damaging  of  several  NWE  environments  affecting  an 
MECHANISM  item 


DOMINANT  KILL  Code 
MECHANISM 
CODE  0 

1 
2 

3 

4 


5 

6 
7 


Meaning 

No  NHE  vulnerability  data  on  item  (or  hard) 
Electromagnetic  pulse  (EMP)  only 
Neutron  fluence  only 
EMP  and  neutron  fluence 

Peak  static  overpressure  (AP)  x  dynamic  pressure 
impulse  (I_)  -  threshold  for  vehicle  overturn 
(K)  only 

APIq  -  K  and  EMP 

APIq  -  K  and  neutron  fluence 

APIg  -  K  and  neutron  fluence  and  EMP 


2.  NUDACC  METHODOLOGY 

2.1  Subroutine  Flow  Chart 


The  current  version  of  IRJDACC  can  perform  a  static  (snapshot) 
evaluation  of  NHE  equipment  damage  and  personnel  casualties  after  a 
nuclear  burst  on  the  battlefield.  All  units  are  considered  stationary; 
the  weapons  are  considered  in  the  order  in  which  their  parameters  are 
entered  as  input  data.  A  narrative  of  the  control  logic  written  in 
structured  form  and  a  subroutine  flow  chart  (fig.  1)  follow. 

•  Read  in  the  unit  location,  size,  and  orientation. 

•  Read  in  the  weapon  location  and  yield. 

•  Calculate  the  maximum  effects  radius. 
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•  Read  in  the  equipment  for  each  unit 


•  Select  emd  detonate  each  weapon* 

For  each  unit, 

•  Calciilate  the  distance  from  the  weapon  burst  point  projected  on  the 

ground  to  the  unit. 

•  If  this  distance  is  less  than  the  maximum  effects  radius  of  the 

weapon,  process  the  unit. 

•  Divide  the  unit  into  grid  squares  50  m  on  a  side,  and  calculate  the 

distance  from  the  weapon  to  the  center  of  each  of  these  grids  (fig. 
2) . 

•  Calculate  the  vaurious  environments  at  the  center  of  each  grid 

square,  and  accumulate  the  dose. 

•  Calculate  the  probability  of  survival  of  all  items  for  which  data 

exist  as  a  result  of  these  environments. 

•  Calculate  the  attrition  of  the  items. 

The  attrition  of  items  (personnel  amd  equipment)  is  calculated 
in  the  following  manner: 

•  If  the  items  have  not  been  distributed  over  the  unit,  then  distribute 

the  items  I  otherwise,  read  the  surviving  items  from  a  random  file. 

•  Multiply  the  items  in  each  grid  square  by  the  appropriate  probability 

of  survival  to  determine  the  number  surviving  at  that  grid. 

•  Sum  each  item  over  the  unit  to  determine  the  number  of  that  particular 

item  surviving  after  that  particular  weapon,  euid  write  those  items  to 
a  file  for  further  processing;  then  redistribute  surviving  personnel 
and  equipment  according  to  the  computed  values  for  each  grid  square, 
and  write  the  distributed  equipment  to  the  random  file. 


The  output  from  NUDACC  is  passed  to  a  SORT  routine,  which 
reorganizes  the  data  for  the  NUDPRIMT  program.  The  output  from  NODPRIMT 
is  a  printout  that  has  em  entry  for  each  unit  and  lists  the  personnel 
and  the  equipment  surviving  each  critical  weapon  and  the  dominant  kill 
oMchanlsm  for  that  item.  The  job  execution  sequence  is  illustrated  in 
figure  3. 
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2.2  Environments 


Sweeney  et  al^  comprehensively  discuss  the  environmental 
eqxiations  used  in  HUDACC.  The  following  environments  are  calculated: 

Radiation  (total  dose) 

Blast  (vehicle  overturn,  -  K) 

Transient  radiation  effects  on  electronics  (TREE)  (neutron 

fluence) 


EMP  (vertical  electric  field) 

The  following  damage  mechanisms  are  considered: 
a.  Personnel  (total  dose) 

( 1 )  Postures  considered 

(a)  Exposed 

(b)  Open  vehicle 

( c )  Foxhole 

(d)  Armored  personnel  carrier 

( e )  Tank 

(2)  Damage  criteria  for  each  posture 

(a)  Immediate  permanent  incapacitation,  undemanding 

tasks  ( 16  krad  [tissue] ) 

(b)  Immediate  permanent  incapacitation,  demanding 

tasks  (8  krad  [tissue] ) 

(c)  Immediate  transient  incapacitaion,  undemanding 

tasks  (3  krad  [tissue]) 

(d)  Immediate  transient  incapacitation,  demanding 

tasks  (2  krad  [tissue]) 

(e)  Latent  lethality  (0.65  krad  [tissue]) 


Sweeney,  Jr,,  C.  Moazed,  and  J.  Wicklund,  Unclear  Weapons 
Environments  for  Vulnerability  Assessments  to  Support  Tactical  nuclear 
Warfare  Studies  (U),  Harry  Diamond  Laboratories  HDL-TM-77-4  (June 
1977),  ( OONFIDEUTIAL) 
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b.  Vehicles  (overturn)— A PI^  -  K 

c.  Electronics 

(1)  Neutron  fluence  (TREE) 

(2)  BMP 

(3)  APIq  -  K  (mounted  in  vehicles) 


Any  items  located  within  the  circle  defined  by  the  projection  of  the 
fireball  on  the  ground  are  considered  destroyed  without  further 
processing. 

Equipment  is  categorized  into  one  of  the  following  nine 
equipment  codes  for  input: 

100-199  VRieeled  vehicles 
200-299  Tracked  vehicles 
300-399  Radios 
400-499  Electronic  equipment 
500-599  Radars  and  sensors 
600-699  Missile  systems 
700-799  Power  generation  equipment 
800-899  Weapon  items 
900-999  Aircraft 

Personnel  and  their  posture  may  also  be  specified  according  to  the 
following  code: 

001  Exposed 

002  Open  vehicle 
003  Foxhole 

004  Armored  personnel  carrier 

005  Tank 

The  number  of  personnel  or  items  of  equipment  is  specified  by  a  three- 
digit  count  immediately  following  the  category  code. 

2.3  Restrictions 

The  current  data  arrays  €md  file  structures  for  both  NUDACC  and 
NUDPRINT  are  limited  to  maxima  of  600  units,  45  equipment  codes  per 
unit,  and  200  weapons. 


3.  INPUT  FILES 


3.1  General 


The  files  described  below,  except  for  the  Equipment  Menu  file, 
are  generated  by  the  user.  The  Equipment  Menu  file  is  directly  related 
to  the  vulnerability  data  organic  to  NUDACC  and  cannot  be  casually 
modified.  All  records  within  the  files  are  expected  to  be  in  card  image 
format . 


3.2  NUDACC  Input  Files 

Unit  data,  data  set  reference  No.  S,  File  1 
Code  FORMAT  ( I7, 3X, 2F10. 0, 3F5. 0 ) 


Card  format:  Column 


Description 


01-07  Unit  identification  number 

08-10  Unused 

11-20  Unit  location  (x  value,  meters) 

21-30  Unit  location  (y  value,  meters) 

31-35  unit  rotation  emgle  (degrees) 

36-40  Unit  length  (meters) 

41-45  Unit  width  (meters) 

46-80  Unused  (or  comments) 


Weapon  data,  data  set  reference  No.  5,  File  2 
Code  FORMAT  (3F10.0) 

Card  format:  Column  Description 

01-10  yield  Ocilotons) 

11-20  Weapon  location  (x  value,  meters) 

21-30  Weapon  location  (y  value,  meters) 

Unit  equipment  file,  data  set  reference  No.  5,  File  3 

Code  FORMAT  (17, 5X,9( 213, IX) ) 


Card  format:  Column 


Description 


01-07  Unit  identification  number 

08-12  Unused 

13-75  9  (equipment  code,  count,  1  space) 

76-80  Unused 


TABU  2.  SORT  ROOTXHB  NASS  8TOKA6B  FIUS 


I 


Data  set 

Routine 

File 

File 

Record  sise 

Moadber  of 

rsfsrsncs  Mo. 

status 

(words) 

records 

SORTIN 

SORT 

s 

OXD 

S 

V 

SORIOOT 

SORT 

s 

HTF 

S 

V 

TABU  3 

.  MUDPRXMT  PROGRAM  NASS  STORAGB 

FIUS 

Data  set 

reference  Mo 

Routine 

• 

File 

type 

File 

statxis 

Record  size 
(words) 

Mmaber  of 
records 

1 

MAIN 

S 

OTD 

1200 

2 

2 

MAIN 

s 

OXD 

5 

V 

11 

MAIM 

R 

OTD 

91 

600 

12 

MAIN 

R 

MID 

12 

200 

5.  JOB  CONTROL  LANGUAGE 

The  run  stream  of  listing  1  executes  NUDACC,  SORT,  and  NUDPRINT  on 
the  HDL  IBM  370/168.  The  FORTRAN  IV  H  Extended  Compiler  is  used. 

LISTING  1.  RUN  STREAM  TO  EXECUTE  NUDACC,  SORT  ROUTINE,  AND  NUDPRINT 
PROGRAM 

//NUDACC  JOB  (Installation  dependent) 

//COMPUTE  EXEC  FORTXCG, FARM. GO- ’MORES' 

//FORT. SYS IN  DD  * 

NUDACC  PROGRAM 

/* 

//GO.FT05F001  DD  * 

UNIT  DATA  DECK 

/* 

//GO.FT05F002  DD  * 

WEAPON  DATA  DECK 

/* 

//GO.FT05F003  DD  * 
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LISTING  1.  RUN  STREAM  TO  EXECUTE  NUDACC,  SORT  ROUTINE,  AND  NUDPRINT 
PROGRAM  (Cont'd) 

UNIT  EQUIPMENT  DECK 

/* 

//GO.FTlOFOOl  DD  DSN-&&RDOSE, DISP-(NEW, DELETE) ,SPACE-(400, (600)) , 

//  DCB-(RECFM-F,LRECL-400) ,UNIT-SYSDA 

//GO.FTllFOOl  DD  DSN-&&RNCDCN, DISP-(NEW, PASS) ,SPACE»(364, (600)) , 

//  DCB-(RECFM-F,LRECL-364) ,UN1T-SYSDA 

//GO.FT12F001  DD  DSN-&&REQPT, DISP-(NEM, DELETE) ,SPACE-(300, (9000)) , 

//  DCB-(RECFM-F,LRECL-300) ,UNIT«SYSDA 

//G0,FT16F001  DD  DSN-&&SORTIN,DISP»(NEW,PASS) ,UNIT-SYSDA, 

//  SPACE-( 1084 , ( 1000 ) ) , DCB-( RECFM-VBS , LRECL-24 , BLKS IZE-1084 ) 

//GO,F17F001  DD  DSN-&&HITCNT,DISP-(NEM,PASS) ,UNIT-V10, 

//  SPACE-(4808,(2)) ,DCB»(RECFM»VS,LRECL»4804) 

//SORT  EXEC  SORT 

//SORTIN  DD  DSN»&&SORTIN,DISP=(OLD, DELETE) 

//SORTOUT  DD  DSN-&&SORTOUT,DISP-(NEW,PASS) ,UNIT-SYSDA 
//  SPACE-(1084 , ( 1000) ) ,DCB»(RECFM-VBS , LRECL-24 , BLKSIZE-1084) 

//SYSIN  DD  * 

SORT  FIELDS-(5.0,4.0,A,13.0,4.0,A,9.0,4.0,A),FORMAT-BI,SIZE-E10000 
END 
/* 

//PRINT  EXEC  FORTXCG 
//FORT. SYSIN  DD  * 

NUDPRINT  PROGRAM 

/* 

//GO.FTOlFOOl  DD  DSN-&&HITCNT,DISP-(OLD, DELETE) 

//GO.FT02F001  DD  DSN-&&SORTOUT,DISP-(OLD, DELETE) 

//GO.FT03F001  DD  * 

UNIT  NAME  DECK 

/* 

//GO.FTlOFOOl  DD  * 


LISTING  1.  RUN  STREAM  TO  EXECUTE  NUDACC,  SORT  ROUTINE,  AND  NUDPRINT 
PROGRAM  (Cont'd) 

EQUIPMENT  NAME  DECK 

/* 

//GO.FTllFOOl  DD  DSN«&&RNCDCN,DISP»(OLD, DELETE) 

//GO.FT12F001  DD  DSN-&&MN, DISP-(NEW, DELETE) ,SPACE-(48, (200)) , 

//  DCB=(RECFM»F,LRECL-48) ,UNIT-VIO 

// 

6.  NUDACC  NAMED  COMMON  BLOCKS 
ATTBLK 

Hie  ATTBLK  block  is  used  to  pass  the  2urray  AT(75,10,10)  between 
subroutines  ATTRIT,  DSTRBT,  RDRER,  SORVIV,  and  WTRER  (fig.  1).  AT  is 
used  as  a  work  area  to  keep  track  of  the  items  of  personnel  emd 
equipment  surviving  over  the  grid  structure  of  the  unit.  After  the 
nupnber  of  items  surviving  a  weapon  has  been  calculated,  the  personnel 
Itons  are  reset  at  their  original  values  due  to  the  cumulative  nature  of 
total  dose. 

DKM 

Ihe  DKM  block  is  used  to  pass  the  array  DKM(3,75)  between 
subroutines  ATTRIT  and  SURVIV.  DKM  is  a  logical  array  used  to  indicate 
the  dominant  kill  mechanism  for  each  equipment  item  over  the  unit. 

DPIQ,  DSBLK,  EMPI,  and  RNF 

The  DPIQ,  DSBLK,  QtPI,  and  RNF  blocks  are  used  to  pass  the  arrays 
DP1Q(10,10),  OOSE(10,10),  a(P1(10,10),  and  RMF(10,10),  respectively, 

between  the  subroutines  ENVCAL  and  PBSURV.  Hie  arrays  respectively 
contain  the  environments  API^,  total  dose,  vertical  electric  field,  and 
neutron  fluence  calculated  an  the  center  of  each  grid  square  of  the  unit 
being  considered. 


NUMHIT 

Hie  NUMHIT  block  is  used  to  pass  the  arrays  NH(600)  and  IC(600) 
betmen  subroutines  CKUNIT,  GOTONK,  and  RDEQPT.  NH  contains  the  number 
of  times  that  a  unit  has  been  within  range  of  a  weapon;  IC  contains  the 
total  number  of  items  associated  with  a  unit.  NH  reflects  the  fact  that 
each  personnel  item  generates  five  survival  numbers.  Hie  NUMHIT  block 
is  eventually  passed  via  temporary  file  to  the  NUDPRINT  program. 


OUTBLK 


The  OUTBLK  block  is  used  to  ^u^range  the  value  NEQP  and  the  arrays 
ICOOE(45)  and  ICOUNT(45)  In  contiguous  areas  of  memory  before  writing  to 
a  random  file.  NEQP  is  the  number  of  equipment  codes  associated  with 
the  unit;  ICODE  contains  the  respective  equipment  codes,  and  ICOUNT 
contains  the  number  of  each  item  of  equipment  represented  by  the  code. 

PSDOSE,  PSDPIQ,  PSEMP1,  and  PSNF 

The  PSDOSE,  PSDPIQ,  PSEMP1,  and  PSNF  blocks  zure  used  to  pass  the 
arrays  PSDOSE(5,4, 10, 10) ,  PSDPIQ( 15, 10, 10) ,  PSEMP1 (2, 10, 10) ,  and 
PSNF(31 ,3, 10, 10) ,  respectively,  between  subroutines  ATTRIT  and  PBSURV. 
The  arrays  respectively  contain  the  probabilities  of  survival  of 
personnel  to  dose,  vehicles  to  airblast  overturn  (API^  ~  K) ,  electronics 
to  vertical  electric  field,  and  electronics  to  neutran  fluence  for  all 
possible  equipment  items  in  each  grid  square  of  a  unit. 

UNTANG 

The  UtfTANG  block  is  used  to  pass  the  array  1HETA(600)  between 
subroutines  DISSQR  and  RDUNIT.  THETA  contains  the  angle  of  rotation  of 
each  unit  about  its  right  rear  corner. 

UNTDIS 

The  UNTDIS  block  is  used  to  pass  the  array  DS(10,10)  between 
subroutines  DISSQR  and  ENVCAL.  DS  contains  the  distance  from  the  weapon 
AGZ  to  the  center  of  each  grid  square  of  the  unit. 

UNTEPT 

The  UNTEPT  block  is  used  to  pass  the  array  IPOINT(45,600)  between 
subroutines  ATTRIT  and  POINT.  IPOINT  contains  the  index  value  of  each 
item  of  nonpersonnel  equipment  in  the  unit.  Ibis  index  is  used  to 
locate  data  relevant  to  the  item  in  arrays  IVLARY,  IDPVUL,  lEPVUL, 
NFIVUL,  and  NFSVUL  found  in  subroutine  ATTRIT. 

UNTFBL 

The  UNTFBL  block  is  used  to  pass  the  array  IFB(10,10)  between 
subroutines  ATTRIT,  ENVCAL,  and  PBSURV.  IFB  is  a  logical  array  that 
indicates  whether  a  grid  square  is  within  the  projection  of  the  fireball 
on  the  ground  and  therefore  destroyed. 
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UNTGRD 


The  ONTGRO  block  Is  used  to  pass  the  eurrays  IL(600)  and  IH(600) 
between  subroutines  GOTONB  and  RDUNIT.  The  arrays  respectively  contain 
the  length  and  the  width  of  the  units  in  number  of  grid  squares. 

UNTHIT 

The  UNTHIT  block  is  used  to  pass  the  array  IH(600)  between 
subroutines  ENVCAL,  GOTONE,  and  RDONIT*  A  nonzero  entry  in  IH  indicates 
that  the  unit  has  been  hit  previously  and  points  to  the  proper  random 
record  containing  accumulated  dose. 

ONTIDN 

The  UNTIDN  block  is  used  to  pass  the  array  IDEN(600)  between 
subroutines  RDEQPT  and  ROUNIT.  IDEN  contains  the  respective  unit 
identification  numbers.  nie  UNTIDN  block  is  passed  via  temporary  file 
to  the  NUDPRINT  program. 

UNTLAW 


The  UNTLAW  block  is  used  to  pass  the  variables  ILEN  and  IWID  between 
subroutines  ATTRIT,  DISSQR,  DSTRBT,  ENVCAL,  GOTONE,  PBSURV,  RDRER, 
SURVIV,  and  WTRER.  These  items  respectively  contain  the  length  amd  the 
width  of  the  unit  under  consideration  in  number  of  grid  squares. 

UNTLOC 

The  UNTLCX:  block  is  used  to  pass  the  arrays  XU(600)  and  YU(600) 
between  subroutines  CKUNIT,  DISSQR,  and  RDUNIT.  The  arrays  respectively 
contain  the  x  2md  y  coordinates  of  the  units  in  meters. 


UNTNCT 

The  UNTNCT  block  is  used  to  pass  the  array  NC(600)  between 
subroutines  DSTRBT,  RDRER,  RDUNIT,  SURVIV,  and  WTRER.  NC  contains  the 
number  of  counts  associated  with  each  unit.  The  Individual  entries 
reflect  the  fact  that  five  counts  are  generated  for  each  personnel  item. 

UNTNDX 

The  UNTNDX  block  is  used  to  pass  the  variables  IRAN  and  IREC  between 
subroutines  ATTRIT,  ENVCAL,  GOTONE,  and  RDUNIT.  IRAN  points  to  the  next 
sequential  random  record  to  which  dose  data  will  be  written  for  a  unit 
exposed  to  nuclear  radiation  for  the  first  time.  IREC  points  to  the 
random  record  containing  dose  information  for  a  unit  that  has  been 
exposed  previously. 
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UNXPER 


The  UMTPER  block  is  used  to  pass  the  array  NPER(600}  between 
subroutines  ATTRIT  and  DSTRBT.  NPBI^  is  the  total  number  of  items  of 
personnel  information  for  a  particular  unit* 

UNTRER 

The  UNTRER  block  is  used  to  pass  the  array  ISTRT(600)  between 
subroutines  RDRER,  ROUNIT,  and  WTRER.  ISTRT  points  to  the  start  of  the 
random  attrited  equipment  records  for  a  particular  unit. 

WEPAGZ 

The  WEPAGZ  block  is  used  to  pass  the  arrays  XA(  200 )  and  yA(  200 ) 
•between  subroutines  CKUNZT  and  RDWEAP.  XA  and  YA  are  the  coordinates  of 
the  AGZ  for  a  particular  weapon. 

WEPNUM 

The  WEPNUM  block  is  used  to  pass  the  variad>le  IWNM  between 
subroutines  CKUNIT,  ENVCAL,  and  SURVZV.  IWNM  is  the  sequence  number  of 
the  weapon  whose  effects  are  being  considered. 

WEPRAD 

The  WEPRAO  block  is  used  to  pass  the  arrays  RN(200),  BTD(200), 

RENP(200),  and  ROPIQ(200)  between  subroutines  LARRAD  and  RDWEAP.  The 
arrays  respectively  contain  the  radius  beyond  which  the  effects  of 
neutron  fluence,  total  dose,  ENP,  and  Are  insignificant  for  a 

particular  weapon. 

WEPYLD 


The  WEPYLD  block  is  used  to  pass  the  array  YIELD(200)  between 
subroutines  CKUNIT  and  RDWEAP.  YIELD  is  the  yield  of  a  particular 
treapon. 


7.  NUDACC  PROGRAM  MODULES 
7.1  Introduction 

NUDACC  was  designed  by  using  the  concept  of  functional 
modules.  Bach  module  resolves  one  of  the  several  problems  facing  an 
analyst  who  wishes  to  know  the  survivability  of  personnel  and  equipment 
items  after  a  weapon  laydown.  For  example,  the  analyst  needs  to  know 


the  location  of  units  (personnel  and  equipment)  and  the  yield  and  the 
AGZ  of  the  various  weapons  f  to  examine  the  magnitude  of  the  several  NWE 
environments  produced  by  each  weapon  at  each  unit  location.  If  the 
magnitude  of  any  one  of  the  environments  is  sufficient  to  damage  the 
unit,  he  calculates  the  probability  of  survival  of  the  various  items 
within  the  unit  due  to  the  environment  and  the  number  of  items  surviving 
the  weapon.  He  also  records  the  total  dose  value  to  idilch  any  future 
dose  from  subsequent  weapon  burst  will  be  added  before  calculating 
subsequent  personnel  survival.  The  MUDACC  program  modules  perform  these 
functions  in  greater  detail. 

7.2  MAIN  Program 

The  MAIN  program  consists  of  four  CALL  statements  to 
subroutines  ROUNIT,  ROWEAP,  ROBQPT,  and  CKUNIT.  The  first  three  of 
these  define  the  problem,  and  the  fourth  performs  the  2malysls. 

RDUNIT 

The  RDUNIT  subroutine  reads  the  location  imd  the  size  of  the 
units  considered  in  the  analysis.  The  number  of  50-m  grid  squares  in 
each  unit  is  calculated,  emd  several  unit  related  2urrays  are  set  to 
zero. 


RDWKAP 

The  ROWEAP  subroutine  reads  the  location  and  the  yield  of  the 
weapons  considered.  The  rafuciraum  radius  of  effect  for  eadi  environment 
is  calculated  by  function  RMAX  emd  stored  in  eurrays. 

RMAX 

(See  ROWEAP.) 

RDEQPT 

The  ROEQPT  subroutine  reads  the  items  of  personnel  and 
equipment  associated  with  a  unit  and  generates  a  random  access  file  of 
equipment  codes  and  counts  for  each  unit. 

CKUNIT 

The  CKUNIT  subroutine  determines  the  largest  of  the  several 
maximum  radii  of  effects  for  a  weapon  with  a  call  to  subroutine  LARRAO 
and  examines  each  unit  to  see  if  it  is  in  range.  A  unit  in  range 
results  in  a  call  to  subroutine  GOTONB. 


7.3  Subroutine  GOTONE 


The  GOTONE  subroutine  locates  the  various  records  required  to 
determine  the  probability  o£  survival  of  the  unit  and  performs  the 
analysis  with  calls  to  subroutines  DISSQR,  ENVCAL,  PBSURV,  and  ATTRIT. 


DISSQR 

The  DISSQR  subroutine  divides  the  unit  into  50-m  squares  and 
computes  the  distance  from  the  weapon  A6Z  to  the  center  of  each  square. 

ENVCAL 


The  ENVCAL  siibroutine  calculates  the  environment  at  the  center 
of  each  50-m  square  after  first  determining  that  the  square  is  outside 
the  fireball  radius  of  the  weapon.  If  the  square  is  inside  the 
fireball,  a  status  switch  is  set  for  reference  during  probability-of- 
survlval  calculations. 

PBSURV 

For  the  environments  at  the  centers  of  each  grid  square  in  the 
unit,  the  PBSURV  subroutine  calculates  probabilities  of  survival  for  all 
items  in  the  unit  that  have  vulnerability  numbers.  Function  PSFUNC 
performs  the  actual  calculation.  These  survival  probabilities  provide  a 
menu  from  which  the  proper  probabllity-of-survival  numbers  are  chosen, 
depending  on  the  unit's  equipment  list. 

PSFUNC 

(See  PBSURV.) 

ATTRIT 

The  ATTRIT  subroutine  calculates  the  number  of  items  of 
personnel  and  equipment  surviving  the  effects  of  the  weapon.  The  items 
of  personnel  and  equipment  are  initially  evenly  distributed  over  the 
unit's  grid  squares,  that  is,  the  first  time  that  a  unit  is  within  range 
of  a  weapon  (hit).  This  information  is  eventually  written  in  to  a 
random  access  file  and  read  directly  if  the  unit  is  subsequently  hit. 
For  items  vulnerable  to  more  than  one  environment,  the  separate 
probabilities  of  survival  are  multiplied,  and  the  product  is  taken  to  be 
the  overall  probability  of  survival  of  the  item. 
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DSTRBT 

The  DSTBBT  subroutine  distributes  the  original  complement  of 
the  unit's  personnel  and  equipment  over  the  unit's  grid  squares.  This 
routine  is  ca)J.ed  only  once:  the  first  time  a  unit  is  hit.  Subroutine 
POINT  is  called  from  this  routine. 

POINT 

The  POINT  subroutine  determines  the  index  value  of  each 
equipment  item  in  the  unit.  The  index  value  is  the  relative  position  of 
the  item's  equipment  code  in  the  arrays  LISTA  emd  LISTS.  This  index 
value  is  the  hey  to  the  various  vulnerability  arrays  in  subroutine 
ATTRIT. 


RDRKR 

The  RDRER  subroutine  reads  the  distributed  items  of  personnel 
and  equipment  surviving  before  the  current  weapon  is  considered.  This 
subroutine  is  called  only  if  the  unit  is  hit  more  than  once. 

SORVIV 

The  SURVIV  subroutine  calculates  the  actual  number  of  items  of 
personnel  and  equipment  surviving  the  current  weapon  by  taking  the  sum 
of  each  item  over  the  unit's  grid  squares*  The  dominant  kill  mechanism 
also  is  determined. 

WTRER 


The  WTRER  subroutine  writes  the  \q>dated  random  access  record 
containing  the  nuodier  of  each  item  of  personnel  and  equipment  surviving 
at  each  grid  square  of  the  unit. 


8.  MODIFYING  NUDAOC  VULNERABILITY  ARRAYS 

The  vulnerability  of  an  item  of  equipment  is  expressed  in  eurrays 
LISTA  and  LISTS  in  subroutine  POINT)  arrays  IVLARA  and  IVLARB,  IDPVLA 
and  IDPVLB,  IBPVLA  and  lEPVLB,  NFIVLA  and  NFIVLB,  and  NFSVLA  and  NFSVLB 
in  subroutine  ATTRIT;  and  arrays  DSNU,  DSSIG,  DPHU.  DPSIG,  KAY,  RNFMO, 
RNFSI6,  BP1MU,  and  EP1SIG  in  subroutine  FBSURV. 

The  arrays  in  subroutines  POINT  and  ATTRIT  do  not  contain  actual 
data,  but  pointers  to  data.  The  similar  geometric  configuration  of 
these  arrays  means  that  information  regarding  any  particular  equipment 
item  will  always  be  located  in  the  same  row  and  column  as  the  item's 
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equipment  number  in  array  LISTA  or  LISTB*  For  exaiiq>le/  equipment  item 
373  is  a  radio  mounted  in  a  vehicle.  First#  locate  the  value  373  in  the 
DATA  statements  for  array  LISTA  or  LISTS  in  subroutine  POINT.  This 
gives  a  row-column  value  of  C3  in  LISTS.  To  determine  the  vulnerability 
information  on  this  item#  locate  row-column  C3  in  IVLASS  in  subroutine 
ATTRIT  and  find  the  value  07.  By  the  dominant  kill  mechanism  code 
(sect.  1),  a  value  of  7  indicates  vulnerability  numbers  for  API^  -  K, 
neutron  fluence#  and  BiP.  To  determine  the  vulnerability  numbers  for 
API^  -  K,  locate  row-column  C3  in  the  DATA  statement  for  array  IDPVLB 
and  find  the  value  05.  This  means  that  the  fifth  elements  in  aurrays 
OPHD#  DPSIG#  and  KAY  in  subroutine  PBSUBV  correspond  to  the  meam  (u)# 
standard  deviation  (a)#  and  K  values  for  API^  -  K  for  item  373.  These 
values  are  1.68,  1.96#  and  8.88#  respectively. 

To  determine  the  vulnerability  numbers  for  neutron  fluence#  go  to 
subroutine  ATTRIT#  locate  row-column  C3  in  the  DATA  statement  for  array 
NFIVLB#  and  find  the  value  01.  In  subroutine  PBSURV  locate  the  first 
elements  in  the  DATA  statements  for  arrays  RNFMU  and  RNFSIG  that 
correspond  to  the  p  and  o  values  for  neutron  fluence  for  item  373;  these 
are  27.3  and  0.37#  respectively.  The  second  element  in  array  T2  gives 
0.9  as  the  neutron  transmission  factor  for  the  vehicle.  Also  of 
Interest  is  the  C3  entry  in  the  DATA  statement  for  array  NFSVLB#  which 
holds  the  value  02.  This  number  points  to  the  neutron  transmission 
factor  for  the  vehicle. 

The  vulnerability  numbers  for  EMP  are  located  in  a  similar  manner. 
Find  row-column  C3  in  the  DATA  statement  for  array  lEPVLB#  which  holds 
the  value  02.  The  second  elements  in  the  DATA  statements  for  arrays 
EPINU  and  EP1SIG  in  subroutine  PBSURV  contain  14.7  and  2.94# 
respectively. 

When  vulnerability  number  updates  are  required  for  item  373#  the  new 
p#  o#  and  K  values  will  have  replaced  the  original  entries.  Since  it  is 
possible  for  several  items  to  use  the  seune  vulnerability  numbers#  assure 
that  the  update  applies  to  all  the  items  whose  pointers  reference  the 
modified  values. 


9.  ADDING  NEW  EQUIPMENT  ITEMS  TO  NUDACC  VULNERABILITY  ARRAYS 

Adding  new  equipment  to  NUDACC  involves  providing  values  to  the 
various  vulnerability  arrays.  The  basic  methodology  outlined  in  section 
8  applies  also  to  locating  the  proper  array  elements  to  be  changed. 


Before  NUDACC  can  be  used  to  calculate  the  probability  of  survival 
of  an  item#  it  is  mandatory  that  p  and  a  be  available  for  the  log-normal 
distribution  representing  the  probability  of  damage  of  the  item  (and  for 
some  items  the  overturn  constant  or  the  transmission  factor). 
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nie  first  step  is  to  select  an  appropriate  equipment  code  for  the 
Item.  This  must  be  one  of  the  numbers  in  the  DATA  statements  for  the 
arrays  LISTA  and  LISTS  In  subroutine  POINT.  Carefully  note  the  row- 
column  location  of  this  number  because  It  will  be  the  key  to  locating 
related  elements  In  the  vulnerability  arrays.  Choose  the  proper  code 
from  the  list  given  In  section  1  to  indicate  the  availability  of 
vulnerability  data  for  the  item.  Enter  this  number  in  the  proper  row- 
column  of  the  DATA  statement  for  array  IVLARA  or  IVLARB  in  subroutine 
ATTRIT. 

Turn  to  subroutine  PBSURV  and  check  the  DATA  statements  for  arrays 
DPMU,  DPSIG,  and  KAY  for  vehicle  overturn;  RNFMU  and  RNFSIG  for  neutron 
fluence  (TREE);  and  EP1MU  amd  EP1SIG  for  EMP,  to  be  certain  that  the 
vxilnerablllty  numbers  for  the  item  are  not  already  Included.  The 
nvimbers  can  be  Included  when  items  that  are  ji^yslcally  similar  have 
different  names  and  model  numbers.  If  the  new  vulnerability  numbers  are 
unique,  increase  the  DIMENSION  statements  for  DPMU.  DPSIG.  and  KAY; 
RNFMU  and  RNFSIG;  and  EP1MU  £uid  EPISIG.  as  appropriate.  Add  the  new 
vulnerability  nuunbers  to  the  DATA  statements  for  these  arrays  and  note 
the  element  numbers  (index  values)  of  the  new  entries.  Then  modify  the 
DO  loops,  which  use  the  data  in  these  arrays  by  increasing  the  test 
values  in  the  DO  statements.  For  APIg  -  K  (vehicle  overturn)  the 
statement  is  this: 

DO  120  M=1.(new  value) 

For  neutron  fluence.  the  statement  is  this: 

DO  210  11^1. (new  value) 

For  EHP.  the  statement  is  this: 

DO  310  11=1. (new  value) 

Return  to  subroutine  ATTRIT  and  modify  the  DATA  statements  for 
arrays  IDPVLA  and  IDPVLB  for  vehicle  overturn.  lEPVLA  and  lEPVLB  for 
EMP.  NFIVLA  and  NFIVLB  for  neutron  fluence.  and  NFSVLA  and  NFSVI£  for 
the  neutron  fluence  transmission  factor  by  entering  the  index  value  of 
the  vulnerability  numbers  in  the  proper  row-column  locations. 

Finally,  add  the  new  equipment  code  and  equipment  name  to  the 
EQUIPMENT  NAME  LIST  file.  This  is  an  input  file  to  the  NUDPRINT  program 
with  data  set  reference  Mo.  10.  Be  careful  to  keep  this  file  in 
ascending  numerical  order. 
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